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Abstract 

Coupled Ginzburg-Landau equations appear in a variety of contexts involving in- 
stabilities in oscillatory media. When the relevant unstable mode is of vectorial 
character (a common situation in nonlinear optics), the pair of coupled equations 
has special symmetries and can be written as a vector complex Ginzburg-Landau 
equation. Dynamical properties of localized structures of topological character in 
this vector-field case are considered. Creation and annihilation processes of differ- 
ent kinds of vector defects are described, and some of them interpreted in theo- 
retical terms. A transition between different regimes of spatiotemporal dynamics is 
described. 

Key words: Vector Ginzburg-Landau equation, topological defects, 
spatiotemporal chaos, optical instabilities, light polarization. 



1 Introduction 



Localized structures, objects with some kind of particle-like behaviour, can 
be found in a variety of non-linearly evolving fields. Examples are vortices 
in fluids, superfluids and superconductors, propagating pulses of excitation in 
nerve systems, solitary waves in chemical media, in parametrically driven sur- 
faces, and in granular media, among others [1] . An understanding of complex 
evolving configurations can be sometimes achieved in terms of the interaction 
rules of the particle-like entities. In most cases, however, numerical integration 

1 http : //www. imedea.uib . es/PhysDept/ 
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of the evolution equations is the most powerful tool to investigate the role of 
the localized structures. 

Nonlinear optical cavities have been specially prolific in providing examples 
of localized structures [2]. They appear in the tranverse profile of the field 
and can take the form of vortices, or of bright or dark dissipative spatial 
solitons. It is usually considered that the polarization degree of freedom of the 
electromagnetic field is fixed either by material anisotropies or by experimental 
arrangement. Thus the description of the dynamics is done in terms of a scalar 
field. However in the cases in which the polarization of the light is not fixed, 
the vector nature of the electromagnetic field leads to striking topological 
phenomena [3]. Recently reported examples of localized structures for which 
the polarization plays a fundamental role include points of zero amplitude in 
both components of the field in a periodic elliptically polarized background, 
found in type II OPO [4] and dots of low amplitude in a circularly polarized 
component of the field in an almost circularly polarized background found in 
self-defocusing vectorial Kerr resonators [5]. 

The scalar Complex Ginzburg-Landau (CGL) equation is considered a para- 
digm model for the qualitative description of general non-linear oscillatory me- 
dia (see, for example, [6-8]). The vector complex Ginzburg-Landau (VCGL) 
equation [9-13] plays the same role when the order parameter is of vector 
character, as is the case for the electromagnetic field when the polarization is 
not fixed. For example, the VCGL equation is an appropriate model for laser 
emission from wide-aperture resonators close to the lasing threshold [13] in the 
absence of polarization-selecting cavity features. Different kinds of localized 
structures [9-12,14-16] are present in the dynamic states of the twodimen- 
sional VCGL equation. From a topological point of view, they are defects, 
that is, places where the field state departs from some basic ordered state. 
These objects carry topological properties which endorse them with a charac- 
teristic stability and robustness. 

In this Paper we study the dynamical properties of these defects for sev- 
eral dynamical regimes of the VCGL equation. In addition to reviewing and 
adding details to results previously presented, which address the synchroniza- 
tion properties of spatiotemporally chaotic states [14], the identification of a 
transition from a glass to a gas phase [15], and the formation and annihila- 
tion processes leading to the different types of defects [16], we interpret some 
of these findings in terms of the stability properties of the waves emitted by 
the defect cores, and in terms of a recent perturbative analytic approach [17]. 
Since we are interested mostly in the dynamics of topological defects, we con- 
sider here just the two dimensional VCGL case, for which point defects are 
topologically stable. For VCGL dynamics in one spatial dimension we refer 
the reader to [18,14,19]. 
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Section 2 presents the basic equations and properties. The different topological 
defects are described in Sect. 3, and Sect. 4 is devoted to their dynamical 
properties. A description of a gas-like phase is presented in Sect. 5. In Sect. 6 we 
briefly discuss a situation of phase separation. Finally (Sect. 7) we summarize 
our conclusions. The Appendix addresses the stability of waves emitted by 
vector defects. 



2 The vector Complex Ginzburg-Landau equation 



Extended systems close to a Hopf bifurcation at zero wavenumber are de- 
scribed by the CGL equation : 

d t A = A+(l + ia)V 2 A-(l + i(3)\A\ 2 A. (1) 



For systems where the relevant unstable mode is of vectorial character, Eq. 
(1) is generalized to the VCGL equation, 



d t A = A + (1 + ia) V 2 A - (1 + if3) 



h - 1) 

^A • A*)A + '-{A ■ A)A* 



•(2) 



In the simplest case A has two complex components. This is precisely the 
case of optical systems, for which, A = (A x ,A y ) describe the complex slowly 
varying amplitude of the electric field, where A x and A y are the cartesian 
components. The right and left circularly polarized components (A + ,A-) are 
related to them via the relations A x = (A + + A_)/y/2 and A y = (A + — 
A_)/(iy/2). In terms of the circular components, Eq. (2) reads 

d t A ± = A± + (1 + ia)V 2 A ± - (1 + i(3){\A ± \ 2 + 7 |^tI 2 )^±- (3) 



The coefficient 7, which in general can be complex, gives the coupling between 
the components. For 7 = we recover a pair of uncoupled equations of the 
form (1). We will often refer to Eq. (1) as the scalar CGL, to contrast it with 
Eqs. (2) or (3). For 7 7^ 0, Eq. (3) can be thought as a particularly symmetric 
example (for example group velocity terms are absent) of a pair of coupled 
complex Ginzburg-Landau equations of the kind usually encountered in wave 
competition situations [6]. 

When Eq. (2) is derived for the Hopf bifurcation leading to laser emission[13], 
a is related to the strength of diffraction, and (3 to the nonlinear frequency 
detuning. In this case, the so called Benjamin- Feir stability criterion 1 + af3 > 
is always satisfied, which means that there are always some plane waves 
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which are stable against long wavelength perturbations. In addition, for laser 
systems, 7 turns out to be a real parameter, and 7 > —1. In this Paper, we 
will restrict our study to these ranges of parameters. 

Equation (3) has a family of plane-wave solutions of the form: 

A± = Q ±e - i ( k ±-"-"±*+*o±) . (4) 

Within this family, the simplest solutions are the circularly polarized traveling 
waves, for which one of the components is identically zero: 

Q + = l-k 2 ,Q_ = 0,u + = -(3 + ((3-a)k 2 (5) 



or 

Q + = 0,Q- = l-k 2 ,Lj_ = -[3 + ([3-a)k 2 (6) 

Within the family (4) there are also solutions in which both components co- 
exist. For wavevectors with the same modulus, k + = k_ = k, we have 



Ql = Q- = Q 2 = 1 



1 + 7 

uj + = uj_ =cu = -f3 + (f3 - a)k 2 . (7) 

We can have a linearly polarized traveling wave (k + = k_) or a standing 
wave with periodic linear polarization (k + = — k_). It is also possible to have 
depolarized solutions [13] with k + 7^ 

Q± = (l- 7 + 74 -*£)/(! -T 2 ) 

u± = -akl -ml + lQ%) (8) 

The qualitative behavior for the VCGL equation (3) for I7I < 1 is rather 
different from the one for 7 > 1. In the second case the nonlinear competition 
between A + and A_ tends to favor one of them against the other, so that in 
regions where A + is developed, A_ generally vanishes, and viceversa. There are 
stable circularly polarized solutions of the form (5) or (6) while solutions of the 
form (7) which involve the coexistence of both fields are linearly unstable. On 
the contrary, for I7I < 1, wave coexistence is generally achieved. Circularly 
polarized solutions of the form (5) or (6) are linearly unstable while there 
are stable solutions of the form (7) and (8). The character of the topological 
defects in both cases is different [11,12]. In Sects. 6-5 we will consider the case 
I7I < I. The case 7 > 1 will be discussed in Sect. 6. 
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For 7 real, the point in parameter space a = (3 = represents a special 
case of particularly simple dynamics. The VCGL equation becomes then a 
potential system, or more precisely [20,21], a 'relaxational gradient' dynamical 
system. Dynamics at long times approaches steady states that are identified 
as the minima of a Lyapunov potential. This is the case which is considered 
in greater detail in [11,12]. A more general situation is a — f3. In this case the 
system can be classified as 'relaxational non-gradient' [20,21], and although 
a Lyapunov functional can be still identified, the attractor for the dynamics 
is no longer steady in general. The limiting 'conservative case' a = (3 — > oo, 
also addressed in [11,12], leads to coupled nonlinear Schrodinger (or Gross- 
Pitaevskii) equations. This last limit is of great interest for the description of 
multicomponent Bose condensates [22] or the propagation of vector solitons 
in optical fibers [23]. 



3 Scalar and vector defects for |7| < 1. 

Starting from random initial conditions, even if the condition for existence 
of stable plane waves (1 + a(3 > 0) is satisfied, solutions of (3) usually do 
not evolve to a simple plane wave (4) due to the presence of defects: points 
where the phase of A + or A_ is undefined and the corresponding amplitude 
is equal to zero. Around each defect, a spiral wave develops (except if a — (3) 
that, far from the defect core, approaches a plane wave with a particular wave 
number which has been dynamically selected by the presence of the defect. For 
example, Fig. 1 shows different views of a configuration obtained after time 
evolution [24] starting from random initial conditions, for a = 0.2, (3 — 2, and 
7 = 0.1. The configuration consist of wave domains separated by shock waves. 
There is a defect in both components (the black dot) in the center of each 
domain. Inhomogeneities and defects in just one of the components have been 
expelled away from the defect core with a given group velocity. 

Therefore, two different kinds of defects are readily identified: Vector defects 
are points at which the two components A + and A_ simultaneously vanish 
(and thus the same happens for A x and A y at the same point). Scalar defects 
(mixed defects in the nomenclature of [11,12]) are those at which only one of 
the two fields A + or A_ vanish. In this case the vector A is not zero, and the 
components A x and A y do not necessarily vanish, but they have still a singular 
topological structure to be described below. The following ansatz describes the 
field A close to a scalar or to a vector defect: 

A±(r,9) = R ± {r)e^ ±{rfi) . (9) 
Where (r, 9) are the polar coordinates with the origin of coordinates at the 
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defect core. The amplitude of A + and (or) go to zero at the core of a vector 
(scalar) defect: R + and (or) R_ — > for r — > 0. The phase of the component 
which vanishes at r = is 

0i(r, 6») = n,j6 + ^i(r) - u{t + O i , (10) 



where the subindex % stands for + or — . From expression (10), rij is seen to 
satisfy 

rii = — fvfa-dr, (11) 



where T is a closed path around the defect. This identifies rij as the topological 
charge associated to the singularity in the field A { . It necessarily takes integer 
values. In our numerical simulations we have never found for it a value different 
from 0, +1 or — 1. uji is the rotation frequency of the spiral wave in A t . For 
r — > oo, ip±(r) - ► 

It is possible to make a classification of defects using topological arguments 
[11,12]. Although this is not the way in which it was originally introduced, 
the classification becomes particularly simple in terms of the charges (11). As 
stated before, a vector defect is the one at which both components of the field 
vanish. Thus, necessarily both n + and n_ are nonvanishing. A vector defect is 
of type argument when the charges in the two field components have the same 
signs, i.e., when n + = n_ = 1 or n + = n_ = — 1. If the charges are of opposite 
signs, i.e., when n + = — n_ = 1 or n + = —n_ = —1, the vectorial defect is of 
director type. A scalar defect is the one at which just one component of the 
field, but not the other, vanish. This implies that only one of the two charges 
is different from zero: n + = ±1 and n_ = 0, or n_ = ±1 and n + = (see 
Fig. 2). 

Numerically we observe that for vectorial defects both fields have identical 
modulus near the core: \A + \ = \A_\, thus R+(r) = R-(r) = R(r). Then close 
to a vector defect, Eq. (3) becomes 

d t A ± = A± + (1 + ia)V 2 A ± - (1 + i(3){\ + i)\A\ 2 ± A ± . (12) 



Each component is described by a scalar CGL equation (1) for a rescaled value 
of the common amplitude, A = y/\ + jA±. In the rescaled variable, the shape 
and the size of the vectorial defect is independent of 7. Furthermore, we can use 
the arguments given by Hagan for the scalar CGL equation [25] to determine 
the wave numbers k + and According to this the spiral wavenumber depends 
only on the parameters a and /3, therefore k + = k_ = k H , where fc^(o;, (5) is the 
wavenumber of a spiral in the scalar CGL. Far from the defect core, as r — > 00, 
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ip+{r) = ip-(r) = ip(r) — > fear. The rotation frequencies are uo + = uo~ = u, 
where oo is given by Eq. (7) using k = kn- It turns out [25] that the dependence 
of kn on a and (3 is always through the combination 



a — (3 



(13) 



K = 



1 + a(3 



a parameter that will be important in the following. If k — 0, kn — 0, so that 
when a = (3 (the relaxational case) a homogeneous state is found instead of a 
wave far from the vector defect core. 

For an argument vectorial defect, n + = n_ = n = ±1, and we have, 



where A0 O — <Po+ — <f>o- is the difference of initial phases of A + and A_. 
If A represents an optical field, its polarization is linear, since the cartesian 
components A x and A y oscillate in phase. The angle of polarization is the 
constant £ = A0 o /2. 

For a director vectorial defect, n + = — n_ = n — ±1, we have, 



The polarization in this case is also linear at each point, but the angle of 
polarization, given by £ = n6 + A0 o /2, rotates around the core of the defect 
(from here the name of director defect). 

The defects just presented do indeed appear spontaneously in numerical simu- 
lations of Eq. (2) in appropriate parameter ranges. However, the classification 
by itself does not indicate much about the dynamics. It turns out that, when 
present, the vectorial defects dominate the dynamics, as shown in Fig. 1. They 
generate a large exclusion domain from which any inhomogeneity or scalar de- 
fect is expelled away. The scalar defects accumulate at the domain limits. The 
modulus of A + and A_ remain essentially frozen in time. The phases <j) + and 
0_ display corotating spirals for the argument defects and counterrotating spi- 
rals for the director defects, in agreement with (10) for the appropriate values 
of the charges rij. We also show in Fig. 1 the global phase <p g = <p + + 0_ and 
the relative phase r = <p + — 0_. Argument and director defects are easily 
distinguished in the plot of the global phase: around an argument defect a 
two-armed spiral is formed (in a closed path around the defect (fi g changes in 




(14) 



A x = V2R(r) cos(n0 + A<j) /2)e i[ ^ r) - u > t+ ^+ + *°- )/2] 
A y = V2R(r) sin(n9 + A(f) /2)e i ^- UJt+ ^ 0++ ' t>0 -^ . 



(15) 
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Att), while a target pattern is seen in the domain of a director defect (<fi g does 
not change). The relative phase <p r does not change in a closed path around 
an argument defect while for a director defect it changes in An. The scalar 
defects are identified as the points where the global or relative phase changes 
by 2tt in a closed path around the point. The modulus and phase of A x are also 
plotted. In the domains of argument defects, the modulus of A x is constant 
and the phase is a spiral, in agreement with Eq. (14). In the domains of the 
director defects, \A X \ has a characteristic radial structure that arises from the 
term cos(n# + A<p /2) in Eq. (15). The phase of A x is a target pattern broken 
by a straight line where there is a phase jump of n. This jump is produced by 
the change of sign of the cosine in Eq. (15). 



Figure 3 shows a scalar defect present in the A + component. It has been 
obtained at the same values of a and (3 as before, but 7 = 0.4. At these pa- 
rameter values, vectorial defects do not appear (as discussed later), so that 
the scalar ones are allowed to generate wave domains, as the one shown in 
Fig. 3. At the scalar defect core, where one of the component vanishes, the 
modulus of the other has a local maximum. Thus the object is circularly po- 
larized and imposes some elliptic polarization to its neighborhood. The phase 
of A + is a spiral wave whose wavenumber can not be derived from the Hagan 
approach, since for scalar defects we have that \A + \ 7^ \A_\ (see Fig. 4) and 
no reduction to a single scalar CGL is possible. The wavenumber of the spiral 
wave is smaller than k H . The phase of the nonvanishing component is almost 
constant. More precisely, this phase slowly changes with radial symmetry as 
one moves away from the center, so that the lines of constant phase form a 
target pattern. In a similar set of equations that include convective terms and 
a complex coupling, scalar defects have also been reported [27]. In that case, 
the phase of the nonvanishing field component presents a clearer target pat- 
tern. Phase target patterns in the scalar CGL equation are also found under 
suitable boundary conditions [28,29] or inhomogeneities [30]. The amplitudes 
and phases of the x and y linearly polarized components of the field around 
a scalar defect, also shown in Figs. 3 and 4, present spiral waves. Note that 
only phase waves, and not amplitude spiral waves, appear in the scalar CGL 
equation. 



Figure 5 shows another view of scalar defects at two values of 7. The maximum 
in the modulus of the nonvanishing component, associated to the presence of 
the defect, is seen to grow with the value of 7, if 7 > 0. This is a consequence 
of the coupling term in Eq. (3), which favors anticorrelation if 7 > (but 
positive correlation if 7 < 0). Notice also that the size of the scalar defect-core 
increases with 7. 



8 



4 Defect dynamics for |7| < 1 



We first discuss all the theoretically possible processes of creation and anni- 
hilation of defects. These processes are illustrated in Fig. 6 where we plot in 
the horizontal axis the topological charge of the A_ component, and in the 
vertical axis the one of A + . Open circles correspond to scalar defects, black 
circles to argument defects and gray circles to director defects, as in Fig. 2. 
The vectorial sum of arrows in the scheme implements the rule of topolog- 
ical charge conservation. The possible processes are: i) creation of a vector 
defect by the coalescence of two scalar defects (illustrated in Fig. 6(a) for an 
argument defect). This process is described in subsection 4.1. ii) Splitting of 
a vector defect in two scalar defects, which is in fact the inverse process, (see 
subsection 4.3). iii) Annihilation of a vector defect by collision with a scalar 
one. Fig. 6(b) represents the annihilation of an argument defect (the black 
dot) by collision with a scalar defect in the A + component (the open dot in 
the vertical axis), giving as a result a scalar defect in the A_ component (the 
open dot in the horizontal). Fig. 6(c) shows a similar annihilation process for 
a director defect (subsection 4.2). 



Aranson and Pismen performed recently [17] an analytical study of the interac- 
tions between scalar defects at different field components, aiming at character- 
izing the stability properties of vector defects: if two scalar defects in different 
components attract each other, a vector defect will be formed in absence of 
other processes, whereas vector defects will be unstable if their components 
repel. A first result is that the interaction is long ranged, in contrast to the 
character of the interaction between defects in the same field. The approach 
in [17] is perturbative for small |7|. The results depend on the value of k, the 
parameter fixing the asymptotic value of the wavenumber. According to [17] 
there is a critical value k c ~ 0.52 such that if k > k c two scalar defects in 
different components attract each other for 7 > [31]. For 7 < 0, two scalar 
defects in different components attract each other if k < k c . We find however 
numerically that for 7 < the vector defects have a short life time. The reason 
is that these defects generate a spiral with a small wavenumber kn so that the 
group velocity at which perturbations are expelled away v g = 2(a — f3)ku (see 
Appendix) is also small. It is then very likely that scalar defects in the neigh- 
borhood of the vector one (excluded from the analysis in [17]) approach its 
vector core and annihilate one of its components. In the following we focus our 
analysis of the dynamical properties of vector defects to the case 7 > 0. Cre- 
ation, annihilation and splitting of vector defects are described in the following 
three subsections. The relaxational case (a = f3) is discussed in subsection 4.4. 
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4-1 Creation of vector defects for 7 > 



If 7 = both components of the field behave independently and vector defects 
are thus not formed. For 7—^0 the density of vector defects goes to zero for 
any value of a and (5. But the behaviour for increasing 7 depends on the values 
of a and f3. It is useful to have in mind the different stability regions in a-f3 
parameter space for the spirals in the scalar CGL equation, as described for 
example in [32]. Spirals simply can be absent (for k — 0), or rather be stable, 
convectively unstable (the spirals remain in place and look stable because 
perturbations are effectively ejected away thanks to the group velocity on the 
spiral wave) or absolutely unstable. When the spirals are absolutely unstable 
defect turbulence develops. The extent of these regions is altered here by the 
coupling between the fields, given by the parameter 7. 

If k 7^ 0, for a coupling < 7 < 7^ (where 7^ depends on a and /3), vectorial 
defects are formed at short times starting from a random initial condition. 
Initially there is a high density of scalar defects. The density decreases as 
defects of opposite charges collide and annihilate in pairs in each component 
of the field. At this point, two scalar defects, one in A + and the other in 
A-, may get close to each other, join and form a bound structure. This later 
stage is reached later as 7 becomes smaller. Figure 7 shows the creation of a 
vectorial defect by displaying the sum of the field amplitudes, \A + \ 2 + |A_| 2 
(a = 0.2, = 2, so that k = 1.29 > k c « 0.52, and 7 = 0.1). For these 
values of a and (3 the scalar CGL (1) is in the region where the wavenumber 
kn of the spirals is convectively unstable [32]. In the temporal sequence it can 
be seen that two dark dots, corresponding to the two scalar defects, one in 
each component of the field, get together and finally form the vector structure 
(the black dot). Immediately after the coalescence of both defects, the group 
velocity around the newly formed vector defect produces an exclusion region 
that precludes the approximation of other scalar defects. The region grows 
until the system arrives to a configuration similar to the one displayed in 
Fig. 1. The system has reached a frozen structure of domains separated by 
shock waves. The configuration in Fig. 1 is qualitatively similar to the ones 
encountered for the scalar CGL. Since very slow relaxation an other similarities 
with structural glasses have been pointed out in that case, the terms glassy 
state or glass phase are applied to it. We will use also these terms for the 
vectorial case (which is also seen to evolve in a very slow time scale). The 
peculiarity here is that, whereas in the scalar CGL the difference between 
the central defect and the defects at its domain border seems to arise from a 
spontaneous symmetry breaking, here all vector defects always expel from its 
domain the scalar ones, so that vector defects, if present, are always found at 
the center of wave domains. 

The glassy configurations occur for relatively small 7. As 7 is increased vec- 
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torial defects become unstable, leading to the melting of the glass phase. We 
have identified two mechanisms by which this process occurs, to which the 
next two subsections are devoted: (i) domain instability and (ii) core insta- 
bility. In the case of domain instability the domain around a vector defect 
becomes unstable and develops inhomogeneities. These irregularities lead to 
paths that permit the approach of a scalar defect that finally collides with 
the core and annihilates one of the components of the vector defect. In the 
core instability case the vector defect simply splits in two scalar defects. This 
process is described in Ref. [12] for the real coefficient case, where a greater 
symmetry between director and argument defects seems to be present and 
both kind of defects become unstable for the same value of 7. As we will see 
in Sect. 4.3, this is not the case for complex coefficients. 

Finally, in the region of parameter space a- (3 where spirals are absolutely un- 
stable [32] for the scalar CGL equation, no vectorial defects are formed. In 
this case we have what is called in the scalar case defect turbulence. Increas- 
ing 7 does not change the qualitative behavior of each component, although 
anticorrelation between the modulus of both components increases. 



4-2 Domain instability 

By 'domain instability' we denote the situation in which the wave domain 
around a vector defect becomes unstable, begins to fluctuate, and develops 
inhomogeneities. Shortly afterwards one of the two charges forming a vector 
defect is annihilated by an external scalar defect. As a result a free scalar 
defect is left in the other component of the field. The region of parameter 
space a-j3 where this process is observed corresponds approximately to the 
region where for the scalar CGL equation the phase spirals are convectively 
unstable [32]. As 7 is increased from zero the stability of the spirals is modified: 
At a given value of 7, the group velocity is not strong enough to overwhelm the 
growth of the perturbations, the spirals becoming absolutely unstable. At this 
point the domains around the vectorial defects are uneffective as exclusion 
zones, so that scalar defects previously confined to the domain border can 
approach the vectorial defect core. Although this picture is valid for both 
kinds of vectorial defects, director defects survive for larger 7 than argument 
ones. For the parameter values of Fig. 1, argument defects become unstable 
at 7 ~ 0.3 (see Fig. 8), while director defects remain stable up to 7 ~ 0.35 
(see Fig. 9). For larger 7 only scalar defects are found numerically. 

The different stability range of argument and director defects can be under- 
stood through a linear stability analysis of the vector spirals focusing in its 
plane-wave structure far from the core. This analysis is performed in the Ap- 
pendix. The main point to be recognized is the difference between the phase 
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structure between argument and director vector defects. In Figs. 10 (a) and 
(b) we show constant phase curves of both components for an argument defect 
and a director defect. The wave vectors of the components A + and A_ in a 
point x, k + (x) and k_(x) are perpendicular to the constant-phase curves. As 
discussed before, k + = = ku, where ku is the wave number of a scalar 
spiral. For the argument defects the difference of wave vectors kR = k + — k 
vanishes at any point of the plane while for the director defects kR 7^ 0. Far 
from the defect core, k R decreases for the director defect, however it remains 
different from within the exclusion island surrounding the defect. 



The linear stability analysis of the wave in the two cases turns out to be 
different. It is shown in the Appendix that for argument defects, a wave with 
k + = k^ = kn becomes unstable if ku > K p , where K p (a,(3,~f) is given in 
Eq. (25), while in the case of director defects this happens if ku > K s , where 
K s (a,/3,j) is given in Eq. (22). 



In figures 11 (a) and (b) we plot stability diagrams in the a- (3 plane, with the 
line 1 + af3 = as a reference. In fig. 11 (a) we show the curves K p = ku 
and K s = k H for 7 = 0.3. To the right of the curves, plane waves with 
k = ku are unstable. In fig. 11 (b) we plot the curve K p = ku for different 
values of 7 (7 = 0,0.3,0.6,0.9,0.99). We can see that the curve moves left 
and upwards as 7 is increased; the same happens with the curve K s = ku- 
Let us consider a point (ceo, A)) such that plane waves with k = ku are stable 
for 7 = 0. As 7 is increased the curve K p = k H crosses the point turning the 
vectorial waves with k R = unstable, this destabilizes the domains of the 
argument defects. If 7 is further increased, then the curve K s = ku crosses 
the point (cko,A)) an d waves with kR 7^ become unstable, destabilizing the 
domains of the director defects. It should be said that the kind of stability 
analysis performed in the Appendix considers extended perturbations. Since 
spiral waves have a group velocity, it may happens that these instabilities are 
of convective type, that is, localized perturbations do not destroy the spiral 
wave but they are advected away from the defect core of the spiral with the 
group velocity [32]. Thus, we can still find stably-looking spirals although the 
corresponding plane waves far from the core are unstable. The calculation of 
the limit of absolute instability to localized perturbations is quite involved, 
but the result for the convective instability suggests that corotating spirals 
become absolutely unstable before counterrotating ones. We believe that the 
analysis presented here is an indication of the kind of instabilities that can 
affect the vectorial spirals and the order in which they appear as the coupling 
parameter 7 is increased. 
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4-3 Core instability 



A core instability of a vector defect occurs when it splits in two scalar defects. 
This process is roughly present in the region of parameter space a-/3 where the 
spirals are stable in the scalar case. The splitting of a director defect is shown 
in Fig. 12 for k = 0.54 and 7 = 0.95. The size of the vectorial defect-core is 
much narrower than the size of the core of the two scalar defects that remain 
at the end of the process. This is a consequence of two facts mentioned before, 
namely that the size of the vectorial defects does not depend on 7, whereas 
the size of scalar defects increases. Also in this eminent defects become 

unstable for smaller 7 than director defects. For example, for the parameter 
values of Fig. 12 argument defects already split for 7 = 0.75. 

After the splitting, both scalar defects may form a bound pair that resembles 
a rotating molecule [17]. In Fig. 13 a temporal sequence is plotted that shows 
a pair of scalar defects that rotate one around the other. For the particular 
values of the figure, however, the distance between the scalar defects increases 
and the angular velocity decreases, thus the molecule finally unbinds. Stable 
molecules, however, are found for other parameter values as reported in [17]. 

4-4 The relaxational case: a = f3 

Qualitative features distinguish the case a — (5 ( k — 0) from the other cases: 
First, no spiral wave is formed around vector defects (since, as stated before, 
kfi — if k — 0). Second, the group velocity with which linear perturbations 
are expelled away from the neighborhood of a vector defect is zero (the group 
velocity is v g = 2(a — /3)kn, as shown in the Appendix). Thus, perturbations 
and other defects are rather free to come arbitrarily close to vector defects if 
k = and we expect them to be frequently annihilated by processes such as 
the ones shown in Fig. 6(b) or (c). 

Starting from random initial conditions we do not find numerically the for- 
mation of vectorial defects in Eq. (3) for a = (3 and 7 > 0. In Fig. 14, the 
quantities \A + \ 2 and |A_| 2 are plotted for a = (3 = and 7 = 0.1. It can be 
seen that there is no superposition of two defects in the same point. Starting 
with an initial condition with a slightly perturbed vector defect, it sponta- 
neously splits into two scalar defects for any positive 7, as can be seen in Fig. 
15 for 7 = 0.2. This behavior is in agreement with the predictions of [17]. 

More in general, approaching the line a = (3, we observe numerically that 
director and argument defects of initial configurations such as the one in Fig. 1 
split spontaneously, even for very small values of 7, again in agreement with 
[17]. 
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5 The gas phase 



For 7 high enough, the vectorial defects always disappear following one of the 
two mechanisms described above. The system then presents a faster disordered 
dynamics compared with the glassy phase with vector defects. It is a kind 
of gas-like phase, dominated by the scalar defects, which are conserved in 
number during very long times. The scalar defects move faster than in the 
glassy phase, in a way reminiscent of atoms in a gas, from which we borrow 
the name. Furthermore, the spiral wavelength around scalar defects increases 
with 7, so that well-developed spirals do not fit in the domains for 7 close to 
one. Domains are thus less effective as exclusion zones and defects even more 
mobile. 

To distinguish the presence of these two different phases we can use the joint 
probability distribution of the modulus of the field components: |^4-|). 
Figure 16 shows a grey scale plot of \A_\) for 7 = 0.1 (glassy phase) 

and 7 = 0.95 (gas phase); in both cases a = 0.2 and (5 — 2. The plot is 
obtained taking the values of (|A + |,|A_|) at different space-time points. For 
7 = 0.1 there is a broad maximum around \A + \ ~ \A_\ ~ 0.85 with deviations 
from these values being rather uncorrelated, except for the points lying in the 
line \A + \ = \A-\. This line shows that the absolute values of both components 
take simultaneously the same value between and 1, situation identifying the 
core of vector defects. For 7 = 0.95 the line \ A + \ = \A_ \ has disappeared due to 
the absence of vector defects. Instead, the probability distribution approaches 
the curve given by \A + \ 2 + |A_| 2 = 1, which indicates anticorrelation between 
\A + \ and \A-\. This anticorrelation is the fingerprint of the dominance of the 
scalar defects (See Figs. 4 and 5). Similar qualitative behaviour is observed 
for other values of a and j3 (given that k > k c ) as 7 is increased. 

The behavior in the glass and in the gas phase can be interpreted in terms of 
synchronization and generalized synchronization, respectively, of spatiotem- 
porally chaotic configurations of the two field components [18,14]. In this con- 
text, another interesting quantity that gives information about the transition 
from the glassy to the gas phase is the mutual information between field 
components, that can be computed from the individual and joint probability 
densities. In general, for two random discrete variables X and Y the mutual 
information 

= -!>(*,,) In ^ 
P(x,y) 

gives a measure of the statistical dependence between both variables, the mu- 
tual information being if and only if X and Y are independent [33]. Numer- 
ical results presented in [14] show that there is a minimum of J(|A + |, \A_\) 
(7 ~ 0.3 at a = 0.2 and (5 — 2), where the variables have been adequately 
discretized, and that this value of 7 can be interpreted as the transition point 
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from the glassy to the gas phase. See [14] and [15] for more details about the 
vortex unbinding transition between the glassy and the gas-like phases. 



6 Dynamics for 7 > 1 

For 7 > 1, as discussed in Sec. 2, linearly polarized states become unstable 
with respect to circularly polarized states. In this case the system segregates 
in two phases, one right circularly polarized and the other left circularly po- 
larized. Figure 17 is an example of a configuration that started from a random 
initial condition (7 = 1.1, a = 0.2, (3 = 2). Phase separation is clearly ob- 
served. As the system evolves, the length of the interface walls is reduced, in 
a domain-coarsening process. Despite the tendency of the concave domains to 
shrink, some defects remain at long times. Vector point defects are no longer 
topologically allowed since the stable uniform states are such that the van- 
ishing of one of the two polarizations in large areas, not just at points, is 
favored. However we still have scalar point defects that are now lying on the 
domains of one circular polarization, and are points at which the correspond- 
ing field amplitude goes to zero. Only one topological charge is defined, the 
associated to the component which vanishes just at one point. They are of the 
type described in [11,12] as having a 'repolarized core' structure: In domains 
filled with one of the polarizations, say +, a defect is a place at which \A + \ 
goes to zero, thus producing a phase singularity (seen in Fig. 17) associated 
to a topological charge n + = ±1; in response to this behavior the compo- 
nent which does not carry the topological charge takes nonzero values in the 
defect-core region. Other structures for the core are in principle possible, as 
for example a 'punched core' structure in which the component that does not 
carry the charge remains vanishing in the defect-core region. As discussed in 
[11,12], topological arguments imply that only the 'repolarized core' should 
be observed for 7 slightly above 1, and this is indeed what Fig. 17 shows. The 
differences in the sizes of the defect cores seen in Fig. 17 disappear at longer 
times. 

We stress that the topological defects in Fig. 17 are different from the localized 
structures found in other systems with coexistence of two equivalent homo- 
geneous states, such as the Swift-Hohenberg equation in adequate parameter 
ranges or the parametrically forced CGL equation. In these systems the oscil- 
lating tails of the fronts connecting the two equivalent homogeneous solutions 
are responsible for maintaining the existence of localized structures against 
the tendency of the domains to shrink [34]. These localized structures are not 
topological defects. Here, on the contrary, it is the topological character of the 
objects what prevents them to disappear after shrinking to a point: defects 
can only disappear by combination with others of opposite charge. 



15 



7 Conclusions 



We have numerically explored the conditions for the existence and stability of 
the different kinds of defects in the two dimensional VCGL equation: vector 
defects of argument or director type, and scalar defects. Dynamic phenomena 
associated with their creation and annihilation processes have been described. 
In particular, we have identified two mechanisms that lead to the destruction 
of the vector defects as the coupling parameter between field components 
is increased: domain instability and core instability. One or the other takes 
place depending on parameters a and /3. For the domain instabilities of vector 
defects we have obtained analytic understanding in terms of the convective 
instability of the emitted spiral wave. For the core instability at a relatively 
large threshold value of 7 there is no available theory so far. 

Our numerical results have been compared with the useful analytical studies 
by Aranson and Pismen [17]. They studied core instabilities as a function of 
a and f3 for small 7. Their results explain our result that in the relaxational 
limit (a = (3) there are no vectorial defects for 7 > 0. They also explain 
the existence of vector defects far from the relaxational limit and for small 7. 
Our main interest, however, was on what happens far from the relaxational 
limit as function of 7. In this regime we have found threshold values of 7 for 
the domain and core instabilities mentioned before. These instabilities are not 
accessible to the perturbative theoretical discussion of Aranson and Pismen 
and there is need of further analytical studies. 

Finally, the dominance of different types of defects at different parameter 
values leads to the identification of different dynamical regimes, such as the 
glassy and the gas phase. 
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Appendix 

We analyze in this Appendix the stability of the wave emitted by vector de- 
fects. Far from the defect core we can approximate the spiral waves by plane 
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waves. The stability of plane waves with wave number kjj will give us infor- 
mation about the stability of the spirals. The stability analysis presented here 
is an extension to two dimensions of the analysis made in one dimension in 
[13] and it is valid for |7| < 1. The effect of a perturbation on the plane wave 
solution (4) can be written as 

A± = (Q + a± ) e - i ( k ±- r -^±*+^o±) ( 16 ) 



where |k+| = |k_| = ku for a vector defect. Q and ou are given by Eq. (7), and 
a± are the complex perturbations. Using the variables 



(17) 







' 3?(a+ + a_) \ 


1 Sj 


... 


S(a + + a ) \ 


I r 




9?(a + — a_) 


[ Ti . 




k S(a + — a_) t 



where dt and $5 mean real and imaginary part, respectively, we obtain lin- 
earized equations for s, sj, r and 77. Instabilities appear first in the phase 
modes. They are identified as the linear combinations of the variables that in 
the space- homogeneous case lead to zero eigenvalues. They are 9 = —(3s + sj 
and ip = —f3r+ri. Performing the change of variables (s, r, sj, 77) — > (s, r, 9, ip), 
and eliminating adiabatically the rapidly decaying amplitude variables s and 
r in terms of the phases, we get the following phase equations: 



d t 9 ={a- P)k s • V0 + (1 + af3) V 2 6 



(1 + /3 2 ) 
2(1 - k 2 ) 



(ks . V ) 2 £ + i±^(k R .V) 2 £ 



1-7 



+(a - /3)k R • 



(1 + /9 2 ) 



(k s - V)(k R -V)V, 



(l-P)(l- 7 ) 
d t xjj = {a- f3)k s ■ W + (1 + a/3)V 2 ^ - (1 + ^ )(1 H " : ; 

(1+/5 2 ) 



+(a - /?)k R • V0 - 



2(l-F)(l- 7 ) 
ks-V)(k R . V)0, 



(k s .V) 2 ^ + T --^(k R -V) 2 ^ 



1 + 7 



(l-P)(l- 7 ) 
where k s = k + + k and k R = k + — k . 

As explained in Sect. 4.2, k R = in the domain surrounding argument defects, 
but k R 7^ 0, although small, for the domain around director defects. In both 
cases one can approximate |kg| ~ 2kn- For small k R , the dominant first-order 
derivative terms in (18) indicate that small perturbations on the waves travel 
at a group velocity v g — 2(a — (3)kn- 



;is) 
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We consider perturbations of small wave vector q. The eigenvalues of the linear 
system (18) are, 



— = — 1 — a (J H 



2Q 2 (l- 7 ' 



;k R ■ q) 2 + (k s ■ q) 2 ] ± ^ 



(19) 



where q — |q|, q = q/g is the direction of the perturbation vector, and S 
aq 2 + ibq — c, with, 



2\2 



{ 7 2 [(k R • q) 2 - (k s • q) 2 ] 2 + 4(k R • q) 2 (k s • q) 2 } 



& = 8( a -/?)^^ y (k R .q) 2 ks.q 
c = 4(a-/5) 2 (k R .q) 2 . 



(20) 



If k R 7^ 0, as is the case for director defects, we can approximate v^S* = 
iy/c(l — 7^q) + 0(q 2 ). The real part of the two eigenvalues are, 



&(A±) 



= -l-a(3 + 



2 ( 1 + P 2 ) h 2 2 Mn( , 
-k cos x± + 0(q), 



Q2(i_ 7 2) 



(21) 



where x± are the angles between the perturbation wave vector q and the 
wave vectors k + and k_. Since |k R | is finite the instability may arise for 
perturbation wavenumber q smaller than |k R |. In this case, considering the 
most dangerous longitudinal perturbations, i.e. when x+ — or X- = 0) we 
obtain a critical wave number K s such that if k > K s the plane waves become 
Eckhaus unstable. Taking 9ft(A±) = in (21) we get, 



Ki 



(1 + aP)(l - 7 ) 



l + a/3)(l-7) + 2(l + /5 2 )' 



(22) 



The equations for d t 9 and d t ip (18) are coupled, and the eigenvector corre- 
sponding to the unstable eigenvalue mixes both variables. This means that 
the sum and the difference of the field perturbations a + and a_ , become un- 
stable simultaneously in this case. 

If k R = 0, as appropriate around argument defects, we have that k + = k , 
and x+ = X- = X, and 



1 + cos 
Q 2 (l-7 2 ) , 



(23) 
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We obtain for the real parts of the eigenvalues, 



g(A±) 
q 2 



= -1 - a(3 + 



2(1 + /? 2 )(1± 7 ) 
Q 2 (l-7 2 ) 



k 2 cos 2 x + 0(q). 



(24) 



We consider again longitudinal perturbations corresponding to x — 0- Taking 
the marginal stability condition in (24) for the most dangerous eigenvalue, 
9ft(A+) = 0, we see that in this case plane waves become unstable for k > K p , 
where K p is given by, 



For k R = 0, the equations (18) for d t 9 and d t ip are decoupled. The most 
restrictive eigenvalue, A + , is associated with the variable ip, which is related 
to the difference of the field perturbations. When if) becomes unstable, the 
difference a + — a_ grows; in optical language this is called a polarization 
instability. 



References 

[1] H. Riecke, in Pattern formation in continuous and coupled systems, ed. by M. 
Golubitsky, D. Luss and S. H. Strogatz (Springer, New York, 1999); L. M. 
Pismen, Vortices in nonlinear fields, (Oxford Univ. Press, New York, 1999). 

[2] F.T. Arecchi et ai. Phys. Rev. Lett. 67, 3751 (1991); K. Staliunas, G. Slekys 
and C. O. Weiss, Phys. Rev. Lett. 79, 2658 (1997); C. O. Weiss et ai. Appl. 
Phys. B 68, 151 (1999); N. N. Rosanov, in Progress in Optics. Vol. 35, Ed. 
by E. Wolf (North-Holland, Amsterdam, 1996); M. Brambilla et ai. Phys. Rev. 
Lett. 79, 2042 (1997); D. Michaelis, U. Peschel and F. Lederer, Phys. Rev. A 
56, R3366 (1997); M. Tlidi, P. Mandel and R. Lefever, Phys. Rev. Lett. 73, 640 
(1994); W. Firth and A.J. Scroggie, Phys. Rev. Lett. 76, 1623 (1996); C. Etrich, 
U. Peschel, and F. Lederer, Phys. Rev. Lett. 79, 2454 (1997); K. Staliunas and 
J.V. Sanchez-Morcillo, Phys. Rev. A 57, 1454 (1998); G.L. Oppo, A.J. Scroggie, 
and W.J. Firth, J. Opt. B: Quant. Semiclass. Opt. 1, 133 (1999). 

[3] R. Bhandari, Phys. Rep 281, 1 (1997). 

[4] M. Santagiustina, E. Hernandez-Garcfa, M. San Miguel, A.J. Scroggie, G.- 
L. Oppo, Polarisation Patterns and Vectorial Defects in Type II Optical 
Parametric Oscillators, submitted (2001). 

[5] R. Gallego, M. San Miguel and R. Toral, Phys. Rev. E 61, 2241 (2000). 

[6] M. C. Cross and P. C. Hohenberg, Rev. Mod. Phys 65, 851 (1993). 




(l + a/?)(l- 7 ) 



(25) 



(1 + «/?)(!- 7 ) + 2(1 +/ 3 2 )(1 + | 7 |)' 



19 



[7] T. Bohr, M. Jensen, G. Paladin, and A. Vulpiani, Dynamical Systems Approach 
to Turbulence (Cambridge Univ. Press, Cambridge, 1998). 

[8] I. S. Aranson, L. Kramer, The World of Complex Ginzburg- Landau Equation, 
preprint (2001). 

[9] L. Gil, Phys. Rev. Lett. 70, 162 (1993). 

[10] L. Gil, Int. J. Bif. Chaos 3, 573 (1993). 

[11] L. M. Pismen, Phys. Rev. Lett. 72, 2557, (1994). 

[12] L. M. Pismen, Physica D 73, 244 (1994). 

[13] M. San Miguel, Phys. Rev. Lett. 75, 425 (1995). 

[14] E. Hernandez-Garcfa, M. Hoyuelos, P. Colet, M. San Miguel and R. Montagne, 
Int. J. Bif. and Chaos 9, 2257 (1999). 

[15] M. Hoyuelos, E. Hernandez-Garcfa, P. Colet and M. San Miguel, Comp. Phys. 
Comm. 121-122, 414 (1999). 

[16] E. Hernandez-Garcfa, M. Hoyuelos, P. Colet, M. San Miguel, Phys. Rev. Lett. 
85, 744 (2000). 

[17] I. S. Aranson and L. M. Pismen, Phys. Rev. Lett. 84, 634 (2000). 

[18] A. Amengual, E. Hernandez-Garcfa, R. Montagne and M. San Miguel, Phys. 
Rev. Lett. 78, 4379 (1997). 

[19] R. Montagne and E. Hernandez-Garcfa, Phys. Lett. A 273, 239 (2000). 

[20] R. Montagne, E. Hernandez-Garcfa, and M. San Miguel, Physica D 96, 47 
(1996). 

[21] M. San Miguel, R. Montagne, A. Amengual, and E. Hernandez-Garcfa, in 
Instabilities and non-equilibrium structutures V, edited by E. Tirapegui and 
W. Zeller (Kluwer, Dordrecht, 1996). 

[22] A.J. Leggett, Rev. Mod. Phys. 73, 307 (2001). 

[23] J. Yang, Y. Tan, Phys. Rev. Lett. 85, 3624 (2000), and references therein. 

[24] The integration algorithm is the one presented in [15], which is a generalization 
to the vector case of the pseudospectral in space and second-order in time 
algorithm for the CGL equation of [26]. We discretize with a 128 x 128 square 
lattice a domain of 64 x 64 space units with periodic boundary conditions. 
Animations of the different dynamic phenomena can be downloaded from 
http : //www. imedea.uib . es/PhysDept/Nonlinear/research_topics/Vcgl2/ 

[25] P. S. Hagan, SIAM J. Appl. Math. 42, 762 (1982). 

[26] R. Montagne, E. Hernandez-Garcfa, A. Amengual, and M. San Miguel, Phys. 
Rev. E 56, 151 (1997). 



20 



[27] J. Lega, C. R. Acad. Sci. Paris 309(11), 1401 (1989). 

[28] V. M. Eguiluz, E. Hernandez-Garcfa and O. Piro, Int. J. of Bif. and Chaos 9, 
2209 (1999). 

[29] V. M. Eguiluz, E. Hernandez-Garcfa and O. Piro, Phys. Rev. E 64, 036205 
(2001). 

[30] M. Hendrey, K. Nam, P. Guzdar, and E. Ott, Phys. Rev. E 62, 7627 (2000). 

[31] Earlier theoretical analysis [11] left open the possibility for the existence of 
vector defects below a critical value of the coupling, 7 C ~ 0.492, for a = /? = 0, 
but the more recent work [17] rules out this possibility, or equivalently 7 C = 0. 

[32] I. S. Aranson, L. Aranson, L. Kramer and A. Weber, Phys. Rev. A 46, 
R2992 (1992); I. S. Aranson, L. Kramer and A. Weber, in Instabilities and 
Nonequilibrium Strutctures IV, pag. 259, edited by E. Tirapegui and W. Zeller 
(Kluwer Academic Publishers, 1993). 

[33] R. J. McEliece, The theory of Information and Coding: a Mathematical 
Framework for Communication, Encyclopedia of Mathematics and its 
Applications, Vol. Ill (Addison- Wesley, New York, 1977). 

[34] D. Gomila, P. Colet, G.-L. Oppo, M. San Miguel, Stable droplets and growth 
laws close to the modulational instability of a domain wall, Phys. Rev. Lett., 
to appear (2001). 



21 




Fig. 1. Vectorial defects for 7 = 0.1, a = 0.2 and (3 = 2 (k = 1.29). From left 
to right, top row: |A + | 2 , |A_| 2 , global phase c/> 9 , and |Ar| 2 - Bottom row: (p + , (ft-, 
relative phase <ft r , and phase of A x . Here and in the following figures, the values are 
coded in gray scale. In the images displaying amplitude values, black means zero 
amplitude, and thus the presence of a defect. 
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Fig. 2. Scheme showing the possible defects according to different combination of the 
charges n + and n_. White dots correspond to scalar defects, black dots to vectorial 
argument defects, and gray dots to vectorial director defects. 
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Fig. 3. Different views of a scalar defect for 7 = 0.4, a = 0.2 and /3 = 2 (k = 1.29), 
present in the component. From left to right, top row: |A+| 2 , |A_| 2 , \A X \ 2 and 
\A y \ 2 ; bottom row: (f) + , 0_, and phases of A x and A y , 
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Fig. 4. Cross sections of the intensities of a scalar defect (same parameters as in Fig. 
3) . The top figure shows the defect in the + component (solid line) , and associated 
maximum in |^4_| 2 (dashed line). Bottom figure: the spiral waves in the x and y 
amplitudes are out of phase (compare with Fig. 3). 
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Fig. 5. Scalar defects at different values of 7, as a function of the two spatial coor- 
dinates. First row: \A+ 1 2 (left) and |A_| 2 (right) at 7 = 0.1. Second row: Idem for 
7 = 0.4. In both cases k = 1.29, and the singular component is A + . 




Fig. 6. Scheme showing some possible creation and destruction processes of defects 
in the n+-n_ plane. See text for details, (a) Creation of an argument defect, (b) 
annihilation of an argument defect due to the collision with a scalar defect with 
charge n_ = — 1, (c) annihilation of a director defect. 
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Fig. 7. Temporal sequence of |^4+| 2 + |^4-| 2 showing the formation of a vectorial 
defect at short times for 7 = 0.1 and k = 1.29. (a) t = 15, (b) t = 30, (c) t = 45, 
and (d) t = 150. 
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Fig. 8. Annihilation of an argument defect, for 7 = 0.3 and k = 1.29. The top figure 
shows the initial condition generated with 7 = 0.1, the box indicates the domain of 
the argument defect, and is amplified in the images below. First row: |^4+| 2 , second 
row: \A_\ 2 . From left to right: t = t + 30, t = t + 60, t = t + 90, and t = t + 120. 
At the final stage scalar defects approach the defect core annihilating the n + charge. 
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Fig. 9. Annihilation of a director defect, for 7 = 0.35 and n = 1.29. The top figure 
shows the initial condition generated with 7 = 0.1, the box indicates the domain of 
the director defect, and is amplified in the images below. First row:|^4 + | 2 , second 
row: \A_\ 2 . From left to right: t = t + 50, t = i + 100, t = t + 200, and t = t + 300. 
As in the argument defect scalar defect approach at the end the defect core, 

and annihilate one of the components of the director defect (the one in A + ). 
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Fig. 10. Lines of constant phase for (a) an argument defect and (b) a director defect. 
Solid curves correspond to 4> + = and it, and dashed curves to 4>- = and ir. We 
can see that in (a) kR = while in (b) kR / 0. 




Fig. 11. Stability diagrams in the a- (3 plane showing different stability regions (the 
unstable regions are located to the right of the curves) . In both figures the dashed 
curve is 1 + a(3 = 0. In figure (a) the curves K p = kn and K s = ku are plotted for 
7 = 0.3 (see Eqs. (25) and (22)). The wave emitted by argument defects is unstable 
to the right of the K p curve. The wave emitted by director defects is unstable to the 
right of the K s curve. In figure (b) we plot the curve K p = ku for different values 
of 7: 0, 0.6, 0.9 and 0.99 as indicated. 
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Fig. 12. Splitting of a director defect, for a = 0.7, (3 = 2 (re = 0.54) and 7 = 0.95. 
The top figure shows the initial condition generated with a value of 7 = 0.9, the 
square indicates the region amplified in the figures below. In the center of the square 
there is a director defect. First row: |^4+| 2 , second row: |A_| 2 . From left to right: 
t = t + 50, t = t + 100, t = t + 150, and t = t Q + 200 
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Fig. 13. Time sequence of a bounded pair of scalar defects after splitting of an argu- 
ment defect, for k = 1.29 and 7 = 0.9. The pair rotates and the bounding distance 
increases with time. The quantity plotted is \A + \ for (a) t = to, (b) t = to + 10, (c) 
t = t + 20, and (b) t = t + 30. 
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Fig. 14. \A + \ 2 and \A-\ 2 for a = 0, (3 = (re = 0, potential limit), 7 = 0.1 and 
t = 100 starting from random initial conditions. There is not spontaneous formation 
of vectorial defects for real coefficients. 




Fig. 15. Temporal sequence showing the splitting of a vectorial defect for real coef- 
ficients (a = (3 = 0) and 7 = 0.2. First row: |^4+| 2 , second row: |^4-| 2 . From left to 
right t = 50, 100, 150 and 200. 
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Fig. 16. Joint probability distribution p(|A^_ |, \A-\) for (a) 7 = 0.1 and (b) 7 = 0.95 
(in both cases k = 1.29). The histograms that are presented in grey levels have 
been constructed by collecting the values of (1^4+1,1^4-1) at all space points of 
the simulation domain, and at 100 temporal samples separated by At = 1. The 
straight diagonal line in the first plot is a signature of the presence of vector de- 
fects (|^4+| = \A-\) whereas the curve in (b) displays the anticorrelation between 
components characteristic of scalar defects. 
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(bottom row) for 7 = 1.1, a = 0.2 and j3 = 2 (k = 1.29), at t = 50. The system 
segregates into two circularly polarized phases, with some defects (the black dots in 
the amplitude plots, associated to the phase singularities in the phase plots). 
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